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Abstract The exact expression of the nuclear structure function describing inclusive 
lepton scattering off nuclei is recalled, and the basic approximations leading to non- 
relativistic and relativistic nuclear y and x scaling, are illustrated. The general and 
systematic features of y-scaling structure functions are pointed out, and a recently 
proposed novel approach to y-scaling, based on a global scaling variable, yc, which 
incorporates the effect of the momentum dependence of the nucleon removal energy, 
and therefore allows the establishment of a direct link between scaling functions and 
momentum distributions, is illustrated and applied to the analysis of a large body of 
data, pertaining to nuclei ranging from the deuteron to Nuclear Matter. A new type 
of scaling phenomenon, the nuclear a;-scaling, based on a proper analysis of nuclear 
quasi-elastic data in terms of the Bjorken scaling variable xb, is shown to occur at 
high values of the four-momentum transfer Q 2 ; the usefulness of nuclear x scaling is 
pointed out. 



1 The nuclear response in inclusive lepton-nucleus 
scattering 

In inclusive lepton scattering off nuclei, A(e, e')X, the nuclear response, or struc- 
ture function, W(v, g 2 ), which represents the deviation of the cross section from 
scattering from free nucleons, has the following exact form: 

/ + OC 7, 2 Z 

fLe^-fc)* < V QiQe-^-^e-^^o > (1) 
-°° 2n i% 

where H is the Hamiltonian of the target nucleus, ^q{Eq) its ground state wave 
function (energy), the momentum operator of nucleon i, and q = |q| and 
v (Q 2 = q 2 — v 2 ) the three-momentum and energy transfers ( for the sake of 
simplicity only the Coulomb interaction has been considered). 

Due to the non commutativity of H and Eq.(l) cannot be evaluated 
exactly. As a matter of fact, one can write: 

= ^ 4-1 + Ei^(U) + (2) 
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and it can be seen that in general: 



£f(u) 

j¥1 



piq 

2M 



^0. 



(3) 



Let us assume, for the time being, that J2j=£i w (l> i)> §m 

replaces the independent kinematical variables q 2 and by g 2 and xo - ■ > w w , 
one obtains from Eq.(l) [1]: 



= 0; if one then 
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and: 
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which shows that at high momentum transfer the reduced function vW(xo,q 2 ) 
should scale to a 5 function. Such a phenomenon will be called non relativistic 
nuclear x-scaling. 

The non relativistic y-scaling is obtained by introducing the quantity [1]: 
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and, assuming again that [if, = ^> one obtains: 

/dk f°° 
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yo) = f(yo) (7) 



where: 



/(y ) - / n(k ± ,k ll )dk ± = 2tt / n(fc)fc 
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(8) 



is the longitudinal momentum distribution and n(k) is the conventional momen- 
tum distribution normalized such that: 



f d 3 kn(k) = f dy f(y ) 
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(9) 



The above picture has to be improved by considering that in reality [if, t^j] ^ 
and that actual experimental data require the use of relativistic kinematics. 
Both improvements will be implemented in the rest of the paper, by means of the 
relativistic plane wave impulse approximation, in which the " final state interac- 
tion" (FSI) term v {^->j) m Eq.(2) is disregarded and the non relativistic 
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quantity (pi+q)/2M is replaced by its relativistic analog \J (pi + q) 2 + M 2 — M. 
One obtains: 



2 



W(is,q 2 ) = J2 <VLi. k Jv|0|^o> • <5(^ + — En — Ea-i) (10) 
/ 

where: kjy = k + q and k = ki are the momenta of the struck nuclcon after 
and before interaction, respectively, En = \/(k + q) 2 + M 2 is the nucleon total 



energy, and E A _i = y M\_ 1 + k 2 the total energy of the final A - 1 system. 

Starting from Eq. (10), we will introduce and discuss the relativistic nuclear y 
and x scaling. 



2 y-scaling 

Eq. (10), shows that the structure function W(y, q 2 ) is governed by the nuclcon 
spectral function P(k, E) which depends on the energy (E) as well as the mo- 
mentum of the nucleons. In the limit of q 2 — > oo it can be shown that [2] : 

Jim qW(is,q 2 ) =F(y) (11) 



OO OO 

F(y) = 2ir J dE J k dkP(k, E) = f(y) — B(y) (12) 

E mi n k min (y,E) 

where: 

oo k min (y,E) 

B(y) = 2TT J dE J kdkP x {k,E) (13) 

E mm \ v \ 

is the so called binding correction, with Pi being that part of P(k, E) generated 
by ground state correlations and f(y) is given by Eq. 8 with y replaced by y. 
The scaling variable y is obtained (see below)from the relativistic energy con- 
servation appearing in Eq. (10) , and represents the longitudinal momentum of 
the nucleons having the minimal value of the removal energy [2] . The interesting 
quantity is f(y), since its knowledge would provide n(k) by inversion of Eq.(8). 
Unfortunately, the extraction of f(y) from the experimental data requires the 
knowledge of both the experimental asymptotic scaling function F(y), and the 
theoretical binding correction B(y). 

Over the past several years there have been vigorous theoretical and exper- 
imental efforts to explore y-scaling over a wide range of nuclei [3], using the 
relativistic scaling variable y, which, recently, has been shown to lead to scaling 
in the relativistic deuteron,within both the light-front [4] and the Bethc-Salpeter 
[5] approaches. 
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In Ref. [2] the asymptotic scaling function F(y) has been obtained by an 
extrapolation procedure of existing data, so as to get rid of FSI, and the longi- 
tudinal momentum distribution f(y) has been obtained by adding to F(y) the 
binding correction B{y) evaluated theoretically. In Ref. [6, 7] it has been as- 
sumed that f(y) obtained in this way represents the experimental longitudinal 
momentum distributions, whose general and universal features, that are essen- 
tially independent of the detailed dynamics, and are valid from A=2 through 
Nuclear Matter, have been pointed out: 

i) /(0) decreases monotonically with A, from ~ lOMeV -1 when A — 2 to 
~ ZMeV' 1 for heavy nuclei; moreover, for y ~ 0, f(y) ~ (a 2 + y 2 )~\ with 
a ranging from ~ AbMeV for A = 2, to ~ 140MeV for A = 56. 

ii) For 50MeV < \y\ < 200MeV, F(y) ~ e~ a2y2 with a ranging from ~ 50MeV 
for A = 2, to ~ WOMeV for A = 56. 

hi) For \y\ > AOOMeV, f(y) ~ C 2 e^ v \ with B ranging from 2.5 x lO^MeV" 1 
for A = 2, to 1 x 10 _3 Mey _1 for A = 56, and, most intriguingly, b = 
6 x 10~ 3 MeV~ 1 , independent of A. 

The following simple form for f(y) yields an excellent representation of these 
general features for all nuclei: 

2 2 

f(y) = x t~l a + G2e ~ m = f° + h ^ 

ot ~~\~ y 

The first term (/ ) dominates the small y-behavior, whereas the second term 
(/i) dominates large y. The systematics of the first term are determined by 
the small and intermediate momentum behaviour of the single particle wave 
function. For \y\ < a this can be straightforwardly understood in terms of a zero 
range approximation and is, therefore, insensitive to details of the microscopic 
dynamics, or of a specific model. The small k behavior of the single particle wave 
function is controlled by its separation energy, (Q = M + Ma-i — Ma = E min ) 
and is given by (k 2 + a 2 ) -1 so a = (2/i<3) 5 , /x being the reduced mass of the 
nucleon. This agrees quantitatively very well with fits to the data summarized 
in (i) above. 

The most intriguing phenomenological characteristic of the data is that f(y) 
falls off exponentially at large y with a similar slope parameter for all nuclei, 
including the deuteron. Since (i) b is almost the same for all nuclei including 
A = 2, i.e., f(y), at large y, appears to be simply the rescaled scaling function 
of the deuteron; and (ii) b(~ 1.18/m) <^ 1/q:_d(« 4.35/m), it can be concluded 
that the term C2e~ 6 ' y ' is related to the short range part of the deuteron wave 
function and reflects the universal nature of NN correlations in nuclei. 

The remaining parameters, C\ and a, can be related to /(0) and the normal- 
ization condition, Eq. (9). Once this is done, there are no adjustable parameters 
for different nuclei. The intermediate range is clearly sensitive to a, the gaussian 
form being dictated by the shell model harmonic oscillator potential. Notice, 
however, that here the gaussian is modulated by the correct \y\ < a behaviour, 
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namely (y 2 + a 2 )~ 1 , thereby ensuring the correct asymptotic wave function. The 
set of parameters of Eq. (14) for various nuclei is presented in Tab. 1. As an 
example, Fig. 1 shows the longitudinal scaling functions f(y) for 2 H, 4 He and 
56 Fe extracted from the experimental data [2] compared to Eq. (14). The fit is 
excellent. 




lQ~ b I i I i i i I i i i I i i i L 
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Figurel. The longitudinal momentum distribution f(y) for 2 H (dotted line), 4 He (full 
line) and 56 Fe (dashed line) corresponding to Eq. (14) with parameters given in Tab.l 
. The points represent the "experimental" f(y) obtained in Ref.[2]. 



Tab lei. The parameters of Eq. (14) for various nuclei 





Ci[MeV] 


a[MeV] 


a[MeV- L ] 


C 2 [MeV- L ] 


b[MeV-'\ 


' 2 H 


18 


45 


6.1 • 10~ ;l 


2.5 ■ 10~ 4 


6- 10" a 


'He 


41 


83 


7.1 • 10~ ;i 


3.3- 10~ 4 


6- 10" a 


4 He 


106 


167 


6.8 • 10 - ' 1 


6.5 ■ 10~ 4 


6- 10 _J 


L'Aq 


83 


166 


5.1 • 10~ ;l 


5.7- 10" 4 


6- 10" a 


bb Fe 


58 


138 


4.6 ■ 10 -3 


6.2 • 10~ 4 


6 ■ 10"'' 



With these observations it is now possible to understand the normaliza- 
tion and evolution of f(y) with A. First note that Eq. (8) implies /(0) = 
if d 3 k nm = (i/ 2 fc) and so is mainly sensitive to small momenta. Since typical 
mean momenta vary from around 50 MeV for the deuteron up to almost 300 MeV 
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for nuclear matter, it is clear why /(0) varies from around 10 for the deuteron 
to around 2-3 for heavy nuclei. More specifically, since C2 <C Ci/a 2 and fx falls 
off so rapidly with y, the normalization integral, Eq. (9), is dominated by small 
y, i.e., by / . This leads to /(0) ~ (7r 1//2 a) _1 = {2tt^Q)^ 1 / 2 which gives an 
excellent fit to the A-dependence of /(0). Since f(y) is constrained by the sum 
rule, Eq. (9), whose normalization is independent of the nucleus, a decrease in 
/(0) as one changes the nucleus must be compensated for by a spreading of the 
curve for larger values of y. Thus, an understanding of f(y) for small y coupled 
with an approximately universal fall-off for large y, together with the constraint 
of the sum rule, leads to an almost model-independent understanding of the gross 
features of the data for all nuclei. 

To sum up, the "experimental" longitudinal momentum distribution can be 
thought of as the incoherent sum of a mean field shell-model contribution, (fo), 
with the correct model-independent small y-behaviour built in, and a "universal" 
deuteron-like correlation contribution (fi). Thus, the momentum distribution, 
n(k), which is obtained from (8), is also a sum of two contributions: n = n + rii. 
This allows a comparison with results from many body calculations in which no 
and ni have been separately calculated. Of particular relevance are not only the 
shapes of no an d fii, but also their normalizations, SWi) = / n o(i)d?k = J fo(i)dy 
which, theoretically, turn out to be, for 4 He, So ~ 0.8 and Si ~ 0.2 [8] whereas 
Eq. (14) yields S = 0.76 and Si = 0.24. 

In order to minimize theoretical uncertainties arising in the subtraction of 
the binding correction B(y) a new scaling variable has been introduced in Ref. 
[6] wich in principle allows a determination of f{y) free of theoretical contami- 
nations. 

The usual scaling variable y is effectively obtained from energy conservation 

v + M A = [(Mx-i + E A _ X ) 2 + k 2 ] 1 / 2 + [M 2 + (k + q) 2 ] 1 / 2 (15) 

by setting k = \y\, — 1, and, most importantly, the excitation energy, 
E* A _ X = 0; thus, y represents the nucleon longitudinal momentum of a nuclcon 
having the minimum value of the removal energy (E = E min ,E A _ 1 = 0). The 
minimum value of the nucleon momentum when q — > 00, becomes k m i n (y, E) = 
\y — (E - E min )\. Only when E = E min does k min (y,E) = \y\, in which case 
S = and F(y) = f(y). However, the final spectator (A— 1) system can be left in 
all possible excited states, including the continuum, so, in general, E* A _ 1 ^ and 
E > E m i n , so B(y) ^ 0, and F(y) ^ /(?/)• Thus, it is the dependence of k m i n on 
E* A _ 1 that gives rise to the binding effect, i.e. to the relation F(y) ^ f(y)- This 
is an unavoidable defect of the usual approach to scaling; as a matter of fact, the 
longitudinal momentum is very different for weakly bound, shell model nuclcons 
(for which E* A _ 1 <~ — 20MeV) and strongly bound, correlated nucleons (for 
which E* A _ X <~ 50 — 200MeV), so that at large values of \y\ the scaling function is 
not related to the longitunal momentum of those nucleons (the strongly bound, 
correlated ones) whose contributions almost entirely exhaust the behaviour of the 
scaling function. In order to establish a general link between experimental data 
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in different regions of the scaling variable, and longitudinal momentum compo- 
nents, one has to conceive a scaling variable which could equally well represent 
longitudinal momenta of both weakly bound and strongly bound nucleons. One 
can account for this in the following way. The large k and E behaviours of the 
Spectral Function are governed by configurations in which the high momentum 
of a correlated nucleon (1, say) is almost entirely balanced by another nucleon 
(2, say), with the spectator (^4 — 2) system taking only a small fraction of k, 
given by the CM momentum of the pair Kqm [9]. Within such a picture, one 
has [9, 10] 

which shows that the excitation energy of the residual nucleus depends both 
upon k and Kcm', if the latter is set equal to zero, the average excitation en- 
ergy for a given value of k is < _E^_ 1 (fc) >= 4^T3I7- By replacing E* A _ 1 in 
Eq. (15) with 4^1517' tne deuteron-like scaling variable yi introduced in [11] 
(see also [12], where the deuteron-like scaling variable was first introduced) is 
obtained, representing the scaling variable pertaining to a "deuteron" with mass 

~ (2) (2) 

M = 2M — E th , where Ej: h ' = \Ea\ — \Ea-2\- Such a scaling variable, how- 
ever, has the unpleasant feature that the effect of the deuteron- like correlations 
are overestimated at low values of yi and , as a result, the correct, shell-model 
picture provided by the usual variable y is lost. When the CM motion of the 
pair is taken into account, such a defect is cured. If the expectation value of 
Eq. (15) is evaluated with realistic spectral functions, for nuclei ranging from 
3 He to Nuclear Matter, one obtaines ([10, 7]) 

< W >= 4 ~ 7 n ] 7 k 2 + b A - c A ^- (17) 
A-12M 2M v ' 

with bA and ca, resulting from the CM motion of the pair, having values ranging 
from UMeV to ASMeV and 3.41 x 10" 1 to 1.66 x 10 _1 , for 3 He and Nuclear 
Matter, respectively. Placing Eq.(17) in Eq.(15) and subtracting the value of the 
average removal energy < E > to counterbalance the value (17) at low values of 
y, a new scaling variable is obtained which effectively takes into account the k- 
dependence of the excitation energy of the residual A— 1 system, both at low and 
high values of y, unlike the usual scaling variable which completely disregards 
E A -n an d the scaling variable 2/2, which overestimate the effects of deuteron-like 
correlations. In the kinematical region of existing experimental data this, global 
scaling variable yo reads as follows [6] 



VG 



q 2 4^ 2 M 2 -W A 4ll/2 



Wa 2 



(18) 



Here, v A = v + M, M = (2A - 3)M/(A - 1) - £ t (2) - (b A + 2M 2 c A - < E >), 
q = q — cava, and W A = ^a 2 q 2 = M 2 + 2vM — Q 2 . For the deuteron 
E\_ x = 0, so y G - y = I - q/2 + [q 2 /4 - (4is d 2 M 2 - Wd^/Wd 2 ] 1 / 2 ) with 
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Vd = v + M d and = va 2 - q 2 = + 2i/M d - Q 2 . For small values of y G , such 

that (j^zj + — c -a-2m ) E ■ > ' * ne USU& 1 variable is recovered. Thus 
y G interpolates between the correlation and the single particle regions. More 
importantly, however, since k min {q,v,E) ~ \y G \, B(y G ) ~ 0, F(y G ) ~ f{y G ). 
Thus, plotting data in terms of y G allows a direct determination of /(j/g)- One 
would therefore expect from the above analysis, the same behaviour of f(y G ) at 
high values of y G for both the deuteron and complex nuclei, unlike what happens 
with the usual scaling function F(y), and the same shell-model behaviour at low 
values of y as predicted by the usual scaling variable. This is, indeed, the case, 
as exhibited in Fig. (2-5), where the direct link between the scaling function 
F(q,y G ) and the longitudinal momentum distributions is manifest. 
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Figure2. The experimental scaling function of 3 He plotted versus the usual, y, 
(crosses), and the global, yc, (open dots) scaling variables, compared with the scaling 
function of H (full dots). The dashed and full lines are the calculated longitudinal 
momentum distributions of 3 He and 2 H respectively. 



We can summarise our conclusions as follows: 



i) The general universal features of the y-scaling function have been identified 
and interpreted in terms of three contributions: a model-independent zero- 
range contribution, a "universal" 2-nucleon correlation contribution and a 
mean field (shell-model) contribution; 

ii) The shape and evolution of the curve have been understood both quantita- 
tively and qualitatively on general grounds; 

iii) A global scaling variable which incorporates the excitation energy of the 




(A — 1) system generated by correlations has been denned, which allows 
one to obtain the longitudinal momentum distributions directly from the 
experimental data without introducing theoretical corrections. In terms of 
this variable the data strongly support the idea that the large y behaviour 
in all nuclei is essentially nothing but a rescaled version of the deuteron. 
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Figure5. The same as in Fig. 2 but for Nuclear Matter. 

3 cc-scaling 

It can be shown [7] that the transition from the non relativistic xo-scaling, dis- 
cussed in section 1, to the relativistic one, can be achieved by using a minimal 
relativity Hamiltonian, i.e. Tj = \/pf + M 2 — M, in which case the variable 
xo = q 2 /2Mv is replaced by the Bjorken scaling variable xq = Q 2 /2Mv . The 
central point, however, is to understand how non relativistic xo-scaling to a 5 
function is affected by the interaction effects due to [H, t^j] ^ 0, and by the 
use of relativistic kinematics. This is discussed in [7] starting from Eqs. (1) and 
(2). The results can be summarised as follows. Starting from the relation: 

uW{x B ,q 2 ) = %W{x B ,q 2 ) = -F(y, q 2 ) (19) 

and expressing y through xb, the following equation, valid around y ~ (xb — 
1) is found: 

vW{x B ,q 2 )^ - 1 a2 +0(Q~ 2 ) (20) 

(XB 1J + M 2 

where we have used for F(y, q 2 ) the expression f(y) cx -i+ a 'i (cf. Eq. 14). 

It can be seen from Eq. (20) that by plotting the experimental data vs. xb, 
one should not expect a 5 function shape, as in x -scaling, but a Lorentzian 
shape with a width decreasing with increasing Q 2 and converging to a finite 
value a 2 /M 2 = 2E min /M, when Q 2 -> oo. The plots shown in Figs. (6)-(9) 
seem indeed to roughly exhibit a Lorentzian behaviour, but experimental data 
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Figure6. The reduced structure function vW(v,q 2 ) oi A He plotted versus the Bjorken 
variable xb- 

at higher values of Q 2 would be necessary to check the prediction of a saturating 
width. 

4 Summary 

The global y-scaling that we have discussed, by identifying the value of the global, 
Da, scaling variable with the value of the nucleon longitudinal momentum, in- 
dependently of the value of the nucleon removal energy, allows one, unlike all 
previous approaches to y-scaling, to establish a direct relation between the scal- 
ing function and the longitudinal momentum distributions. Nuclear cc-scaling, 
i.e. scaling of inclusive quasi elastic data on nuclei, when plotted vs. the Bjorken 
scaling variable x b , has been shown to qualitatively occur; it can provide useful 
and complementary information on the nucleon momentum distributions. 
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Figure7. The log plot of Fig.6. 
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Figure8. The same as in Fig.6 but for 56 Fe. 
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